Entanglement of Dirac fields in an expanding spacetime 
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We study the entanglement generated between Dirac modes in a 2-dimensional conformally 
flat Robertson- Walker universe. We find radical qualitative differences between the bosonic and 
fermionic entanglement generated by the expansion. The particular way in which fermionic fields 
get entangled encodes more information about the underlying space-time than the bosonic case, 
thereby allowing us to reconstruct the parameters of the history of the expansion. This highlights 
the importance of bosonic/fermionic statistics to account for relativistic effects on the entanglement 
of quantum fields. 
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I. INTRODUCTION 

The phenomenon of entanglement has been extensively 
studied in non-relativistic settings. Much of the interest 
on this quantum property has stemmed from its relevance 
in quantum information theory. However, relatively little 
is known about relativistic effects on entanglement [1-21] 
despite the fact that many of the systems used in the im- 
plementation of quantum information involve relativistic 
systems such as photons. The vast majority of investiga- 
tions on entanglement assume that the world is flat and 
non-relativistic. Understanding entanglement in space- 
time is ultimately necessary because the world is fun- 
damentally relativistic. Moreover, entanglement plays a 
prominent role in black hole thermodynamics [22-29] and 
in the information loss problem [11, 30-34]. 

Recently, there has been increased interest in under- 
standing entanglement and quantum communication in 
black hole spacetimes [21, 35-37] and in using quantum 
information techniques to address questions in gravity 
[38, 39]. Studies on relativistic entanglement show that 
conceptually important qualitative differences to a non- 
relativistic treatment arise. For instance, entanglement 
was found to be an observer-dependent property that is 
degraded from the perspective of accelerated observers 
moving in flat spacetime [6, 7, 9, 16, 40]. These results 
suggest that entanglement in curved spacetime might not 
be an invariant concept. 

In this paper we study the creation of entanglement be- 
tween Dirac modes due to the expansion of a Robertson- 
Walker spacetime. A general study of entanglement in 
curved spacetime is problematic because particle states 
cannot always be defined in a meaningful way. How- 
ever, it has been possible to learn about certain aspects 
of entanglement in curved spacetimes that have asymp- 
totically flat regions [4, 5, 34, 41]. Such studies show 
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that entanglement can be created by the dynamics of the 
underlying spacetime [15, 41] as well as destroyed by the 
loss of information in the presence of a spacetime horizon 
[6, 21, 34]. 

Such investigations not only deepen our understand- 
ing of entanglement but also offer the prospect of em- 
ploying entanglement as a tool to learn about curved 
spacetime. For example, the entanglement generated be- 
tween bosonic modes due to the expansion of a model 
2-dimensional universe was shown to contain information 
about its history [41], affording the possibility of deduc- 
ing cosmological parameters of the underlying spacetime 
from the entanglement. This novel way of obtaining in- 
formation about cosmological parameters could provide 
new insight into the early universe both theoretically 
(incorporating into cosmology entanglement as a purely 
quantum effect produced by gravitational interactions in 
an expanding universe) and experimentally (either by de- 
velopment of methods to measure entanglement between 
modes of the background fields or by measuring entangle- 
ment creation in condensed matter analogs of expanding 
space-time [42, 43]). Other interesting results show that 
entanglement plays a role in the thermodynamic proper- 
ties of Robertson- Walker type spacetimes [44] and can in 
principle be used to distinguish between different space- 
times [15] and probe spacetime fluctuations [45]. 

Here we consider entanglement between modes of a 
Dirac field in a 2-dimensional Robertson- Walker uni- 
verse. We find that the entanglement generated by the 
expansion of the universe for the same fixed conditions 
is lower than for the bosonic case [41]. However we also 
find that fermionic entanglement codifies more informa- 
tion about the underlying spacetime structure. These 
contrasts are commensurate with the fiat spacetime case, 
in which entanglement in fermionic systems was found to 
be more robust against acceleration than that in bosonic 
systems [7, 16]. In the limit of infinite acceleration 
fermionic entanglement remains finite due to statistical 
effects [46, 47] which resemble those found in first quan- 
tization scenarios [48] . 



Our paper is organized as follows: In section II we 
revise the Dirac equation in a spatially flat d-dimensional 
Robertson- Walker universe. Subsequently setting d — 2, 
in section III we calculate the entanglement entropy of 
two Dirac modes and compare it with the bosonic case. 
In section IV we explain the origin of the entanglement 
peculiarities of the fermionic case, showing how it can 
give us more information about the parameters of the 
expansion. Conclusions are presented in section V. 



II. DIRAC FIELD IN A d-DIMENSIONAL 
ROBERTSON- WALKER UNIVERSE 

As we mentioned before, entanglement between modes 
of a quantum field in curved spacetime can be investi- 
gated in special cases where the spacetime has at least 
two asymptotically flat regions. Such is the case of the 
Robertson- Walker universe where spacetime is flat in the 
distant past and in the far future. In this section, fol- 
lowing the work done by Bernard and Duncan [49, 50], 
we find the state of a Dirac field in the far future that 
corresponds to a vacuum state in the remote past. 

Consider a Dirac field ip with mass m on a, d- 
dimensional spatially flat Robertson- Walker spacetime 
with line element, 



ds^ = C{ri){-dr]'^ + dxidx'). 



(1) 



Xi are the spacial coordinates and the temporal coordi- 
nate T] is called the conformal time to distinguish it from 
the cosmological time t. The metric is conformally fiat, as 
are all Robertson- Walker metrics. The dynamics of the 
field is given by the covariant form of the Dirac equation 
on a curved background, 



{ij^(d^-r^) + m}ip = o, 



(2) 



where 7*^ are the curved Dirac-Pauli matrices and F^ 
are spinorial affine connections. The curved Dirac-Pauli 
matrices satisfy the condition. 
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where g^'^ is the spacetime metric. In the flat case where 
the metric is given by 7]"^^ , the constant special relativis- 
tic matrices are defined by. 
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(4) 



The relation between curved and flat 7 matrices is given 
by 7/j — ^iTla where e " is the vierbein (tetrad) field 
satisfying the relation e^ej^rja/j — g^i^- 

In order to find the solutions to the Dirac equation 
Eq. (2) on this spacetime, we note that C{rj) is indepen- 
dent of X. We exploit the resulting spatial translational 
invariance and separate the solutions into 

(5) 



where fc^ — |fep = J2i=i ^f- Inserting this into the Dirac 
equation, we obtain the following coupled equations 



d^ + w?C ± imCC-^l'^ + |fc|' 



€' 



0, 



(6) 



using the fact that the eigenvalues of 7° are ±1. In 
order to quantize the field and express it in terms of 
creation and annihilation operators, positive and nega- 
tive frequency modes must be identified. This cannot be 
done globally. However positive and negative frequency 
modes can be identified in the far past and future where 
the spacetime admits timclikc killing vector fields ^d/drj. 
Provided (7(77) is constant in the far past 77 -^ — cxd and 
far future rj — > +00, the asymptotic solutions of Eq. (6) 
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(7) 



fJ-out = my/C'{+(Xj). 

The action of the Killing vector field on the asymptotic 
solutions allow us to identify (jy^^ * and (jy^l* as negative 
frequency solutions. The sign flip is due to the explicit 
factor i in (6). A consequence of the linear transforma- 
tion properties of such functions is that the Bogolubov 
transformations associated with the transformation be- 
tween in and out solutions take the simple form [50] 
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where aj and (3^ are Bogoliubov coefficients. 

The curved-space spinor solutions of the Dirac equa- 
tion are defined by (with corresponding Uout, Vout and 

Kout), 



C/,„(fc, A;x,r/) = i^„(fc)[C(77)](i-'^)/4 
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where K^^ = -(l/|fc|)((wi„ - /i„)/2^i„)i/2 g^^^ u(0. A), 
w(0. A) are fiat space spinors satisfying, 

7%(0,A) == -m(0,A), 

7%(0,A) = iv{0,X), 
for I < A < 2^^/2-1^ rpj^g ggy ij^ ^j^g "jjj" region can then 
be expanded as. 
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with a similar expression for the "out" region. The in and 
out creation and annihilation operators for particles and 
anti-particles obey the usual anticommutation relations. 
Using the Bogoliubov transformation one can expand the 
out operators in terms of in operators 



aout{k,X) = 



^inf-^out 



K. 



out 



al ain{k,\) 



bout{k,X) = 



^)*^X,,,(-fe)feL(-fe,A') , (11) 
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This yields the following relationship between Bogoliubov 
coefficients, 
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We consider the special solvable case presented in [50] 
C(?7) = (1 + e(l + tanh/9?7))^, where e, p are positive real 
parameters controlling the total volume and rapidity of 
the expansion, respectively. In this case the solutions of 
the Dirac equation that in remote past reduce to positive 
frequency modes are, 
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where Fi is the ordinary hypergeometric function. Sim- 
ilarly, one may find a complete set of modes of the field 
that behaving as positive and negative frequency modes 
in the far future, 
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where uj± — {uJout iwi„)/2. The spacetime obtained by 
considering this special form of C{vi) was introduced by 
Duncan [50]. It is easy to see that it corresponds to a 



Minkowskian spacetime in the far future and past, i.e., 
C — > (1 -I- 2e)^ in the out region and C — t- 1 at the in 
region. 



If we define [7 [^ 
we get that 
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An analogous procedure can be followed for scalar 
fields [41]. The time dependent Klein-Gordon equation 
in this spacetime is given by 



{dl + e + C{^)m')xk{ri)^Q. 



(17) 



After some algebra, the solutions of Klein-Gordon equa- 
tion behaving as positive frequency modes as r; — > 
— oo(i — ?• —00), are found to be 
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Similarly we have 
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where w = {m?{2e+\)'^ — p^)^/"^ . Computing the quotient 
of the Bogoliubov coefficients for this bosonic case, we 
find 



\l-B? 



cosh -Lo 
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cosh -Lu + cosh — CJ+ 



(18) 



III. ENTANGLEMENT GENERATED DUE TO 
THE EXPANSION OF THE UNIVERSE 

It is then possible to find the state in the far future 
that corresponds to the vacuum state in the far past. By 
doing that we will show that the vacuum state of the field 
in the asymptotic past evolves to an entangled state in 
the asymptotic future. The entanglement generated by 
the expansion codifies information about the parameters 
of the expansion, this information is more easily obtained 
from fermionic fields than bosonic, as we will show below. 

Since we want to study fundamental behaviour we will 
consider the 2-dimensional case, which has all the funda- 
mental features of the higher dimensional settings. 



Using the relationship between particle operators in 
asymptotic times, 



brnik) = a^ *hout{k) + ^^ * x{k)al^ti-k) 



(19) 



we can obtain the asymptotically past vacuum state in 
terms of the asymptotically future Fock basis. Demand- 
ing that 6i„(fe,A)|0)m = we can find the "in" vacuum 
in terms of the "out" modes. Due to the form of the 
Bogoliubov transformations the "in" vacuum must be of 
the form 



= J|(Ao|0)o«t + Ai\lkl-k)out) 



where to compress notation |l_fe) represents an antipar- 
ticle mode with momentum ~k and |lfc) a particle mode 
with momentum k. Here we wrote the state for each fre- 
quency in the Schmidt decomposition. Since different k 
do not mix it is enough to consider only one frequency. 
Imposing 6i„(fc)|0)i„ = we obtain the following condi- 
tion on the vacuum coefficients 



.^*Ai|i_0 + /Sfe"*x(fcMo|i-fc) = o 



givmg 



where 



Ai 



K 



^x(A:)Ao = -7-*x(fc)Ao 



r*(k) = 



Pk 



(20) 



(21) 



(22) 



From the vacuum normalization, 

1 = ^n{m^n = \Ao?{l + \l-{k)x{k)?)■ (23) 

Therefore, the vacuum state 

V ^i + b-{k)x{k)? 

is an entangled state of particle modes and antiparticle 
modes with opposite momenta. 

Since the state is pure, the entanglement is quanti- 
fied by the von-Neumann entropy given by S{pk) — 
Tr(pfc log2 Pk) where pk is the reduced density matrix 
of the state for mode k. Tracing over the antiparti- 
cle modes with momentum — fc (or alternatively, particle 
modes with momentum k) we obtain 
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Using the following identity 

XAaI ~ '^P'OutKg^i{p,out ^ ^out) - 

we obtain the entanglement entropy 
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Using (28) we find that the fermionic entanglement is 



/ 



Sf = log 



\ 



I7i.| 






WlF 



(29) 



where \'jp\ — \'^kX{k)\- Note that for massless fields 
(to = 0) the entanglement vanishes since W- — and 
7~=0. Comparing our result to the bosonic case studied 
in [41] we find 
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/ 



where the expression for 73 in (18) differs from that in 
ref. [41] due to the different scale factor used. 

The difference between the bosonic and fermionic cases 
means that the response of entanglement to the dynam- 
ics of the expansion of the universe depends on the na- 
ture of the quantum field. We see from (24) that each 
fermionic field mode is always in a qubit state (the ex- 
clusion principle imposes a dimension-2 Hilbert space for 
the partial state). However in the bosonic case [41] the 
Hilbert space for each mode is of infinite dimension, as 
every occupation number state of the out Fock basis par- 
ticipates in the in vacuum. In both cases the entangle- 
ment increases monotonically with the expansion rate p 
and the total volume expansion parameter e. It is possi- 
ble to find analytically the asymptotic values that both 
fermionic and bosonic entanglement reach at infinity. For 
example, when k = m = p = 1 wc find that as e — ;■ cxd 



7b ^ e--^' 
respectively yielding 



•7r\/2 _ TT 

^p ^ e-^^ — (31) 



5*5(6^00) w 0.0913 S'F(e ^ 00) « 0.0048 (32) 



The entanglement entropy is bounded by Se < log2 N 
where N is the Hilbert space dimension of the partial 
state. The fermionic upper limit Se = ^ corresponds to 
a maximally entangled state. For bosons the unbounded 
dimension of the Hilbert space implies the entropy of 
entanglement is not bounded by unity [41]. However this 
distinction does not guarantee that we can extract more 
information from bosons as we shall now demonstrate. 



IV. FERMIONIC ENTANGLEMENT AND THE 
EXPANSION OF THE UNIVERSE 

As seen in figures 1 and 2 the entanglement behaviour 
is completely different for bosons (fig. 1) and fermions 
(fig. 2). Although the behaviour as the mass of the field 
varies seems qualitatively similar, the variation with the 
frequency of the mode is completely different. 

The entanglement dependence on |fc| for bosons is 
monotonically decreasing whereas for fermions, the 
global space-time structure 'selects' one value of |fc| for 
which the expansion of the space-time generates a larger 
amount of entanglement (peak in figure 2). We shall 
see that this selection of a privileged mode is sensitive 
to the expansion parameters. This may be related to 
the fermionic nature of the field insofar as the exclusion 
principle impedes entanglement for too small |A;|. 

Regardless of its origin, we can take advantage of 
this special behaviour for fermionic fields to use the 
expansion-generated entanglement to engineer a method 
to obtain information about the underlying space-time 
more efficiently than for bosons. 
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Using fermionic fields to extract information 
from the ST structure 



Doing a conjoint analysis of the mass and momentum 
dependence of the entropy we can exploit the character- 
istic peak that SsiiTi, \k\) presents for fermionic fields to 
obtain information from the underlying structure of the 
space-time better than we can do with a bosonic field. 
Let us first show both dependences simultaneously. Fig- 
ure 3 shows the entropy of entanglement as a function of 
our field parameters (|fc| and mass) for different values of 
the rapidity. 

We see from figure 3 that there is no saturation as 
p — >■ cx). Instead, as p is increased the plot is just rescaled. 
This is crucial in order to be able to trace back the met- 
ric parameters from entanglement creation. We also see 
from figure 3 that, for a given field mass, there is an op- 
timal value of |fe| that maximises the entropy. In figure 
4 we represent this optimal \k\ as a function of the mass 
for different values of p, showing how the mode which 
get most entangled as a result of the spacetime expan- 
sion changes with the mass field for different rapidities. 



FIG. 1: Bosonic field: Se for a fixed mass m = 1 as a funciton 
of jfej (up) and for a fixed |fc| = 1 as a function of m (down) 
for different rapidities p = 1, . . . , 100. An asymptotic regime 
is reached when p — > oo. e is fixed e — 1 



From the figure we can readily notice two important 
features 

• The optimal |fc| curve is very sensitive to p varia- 
tions and there is no saturation (no accumulation 
of these lines) as p is increased. 

• There is always a field mass for which the optimal 
|fe| clearly distinguishes arbitrarily large values of 
P- 

In figure 5 we can see a consequence of the re-scaling 
(instead of saturation) of SEijn, \k\) when p varies. In 
this figure we show simultaneously the entropy in the op- 
timal curve and the value of the optimal |fc| as a function 
of the mass of the field for two different values of p, show- 
ing that if p results to be very large, entanglement decays 
more slowly for higher masses. 




FIG. 2: Fermionic field: Se for a fixed mass m = 1 as a 
funciton of |fc| (up) and for a fixed |A;| = 1 as a function of m 
(down) for different rapidities p = 1, . . . , 100. Tlie maximum 
shifts to tlie riglit reacliing an asymptotic value as p — !> oo. e 
is fixed e = 1. The behaviour as |fe| varies is radically different 
from the bosonic case. 




FIG. 3: (Color online) Ssim, \k\) for e = 1, p = 1 (up) p = 
100 (down). Red color = higher SEijn, \k\) 



B. Optimal |fc| tuning method 



1. Part I: Rapidity estimation protocol 



The relationship between mass and the optimal fre- 
quency is very sensitive to variations in p, presenting no 
saturation. Conversely figure 6 shows that the optimal 
A; I curve is almost completely insensitive to e. All the 
different e curves are very close to each other. We can 
take advantage of this to estimate the rapidity indepen- 
dently of the value of e using the entanglement induced 
by the expansion on fermionic fields. 



Given a field of fixed mass, we obtain the entangle- 
ment for different modes fci , . . . , fc„ of the field. Then 
the mode ki that returns the maximum entropy will cod- 
ify information about the rapidity p, as seen in figure 4. 
One advantage of this method is that there is no need to 
assume a fixed e to estimate p, since the tuning curves 
(fig. 4) are have low sensitivity to e (fig. 6). Further- 
more this method does not saturate for higher values of 
p since we can use heavier fields to overcome the satura- 
tion observed in figure 2. While one might expect that 
heavier masses would mean smaller maximum entropy, 
figure 5 shows that if p is high enough to force us to look 
at heavier fields to improve its estimation, the amount of 



bU 








^ 




0.4 


50 




^ 




-- 




0.35 
0.3 


40 




/^^^ 


p=50^^ 


- 




0.25 


^ 30 


/; 


/>-^^ 


9^'ia,,..---' 


- 




c/," 0.2 


20 


£^ 


^-^^'^ 


M£_^_— 


-^ 




0.15 

0.1 


10 


i^ 


-— 


__— t?— - 






0.05 


) 20 


40 60 


80 


100 






mass 








900 












0.4 








..-^ 




800 






p=2000,,-'^^^ 


- 




0.35 


700 




/- 


^^ ^=^Wia,.---^ 


- 




0.3 


600 




X^ 


p=1000,^ 


"^ 




0.25 


500 
400 


/ 
// 


^y^^^ 


p=500^__. 


-- 




kT 0.2 


300 


/> 


"^^^^-^^ 




. 




0.15 


200 






p=100 


. 




0.1 


100 


■ /|S— - 


— ^ — ^ 




- 




0.05 





r"^ 














200 400 600 800 1000 
mass 




400 600 

mass 



FIG. 4: e = 1, optimal |fc| curves (maximum entanglement 
mode) as a function of the field mass for p = 10, . . . , 2000 



FIG. 5: (Color online) Se (blue continuous line) and k (red 
dashed line) in the line of optimal fc as a function of mass for 
p = 10 (up) and p = 1500 (down), e is fixed e = 1 



entanglement will also be high enough due to the scal- 
ing properties of S'b(|A;|, m). We can therefore safely use 
more massive fields to do estimate p since they better 
codify its value. 

Hence we have a method for extracting information 
about p that is not affected by the value of e. Information 
about p is quite clearly encoded in the optimal |fc| curve, 
which is a direct consequence of the peaked behaviour of 
S'_E(|fc|,m). 



2. Part II: Lower bound for e via optimal \k\ tuning 

We can see from figure 4 that for different values of p 
the maximum value for the entanglement at the optimal 
point (optimal k and optimal m) is always S'™'*'^ « 0.35. 
Consider now e ^ 1. In figure 7 we can see how the 
maximum entanglement that can be achieved for optimal 



frequency and mass varies with the volume parameter e. 
Indeed, the maximum possible entanglement that the op- 
timal mode can achieve is a function of only e and is in- 
dependent of p. Hence information about e is encoded in 
the maximum achievable fermionic entanglement. Conse- 
quently we can find a method for obtaining a lower bound 
for the total volume of the expansion of the space-time 
regardless of the value of the rapidity. 

In this fashion we obtain a lower bound for e since the 
entanglement measured for the optimal mode is never 
larger than the maximum achievable entanglement rep- 
resented in figure 7, SE{\k\,m) < S'^^^. For instance 
if the entanglement in the optimal mode is Se > 0.35 
this will tell us that e > 1, whereas if Se > 0.87 then 
we can infer that e > 10. Note that as e increases the 
entanglement in the optimal |fc| mode for the optimal 
mass field approaches that of a maximally entangled state 
when e ^^ c» . 




FIG. 6: (Color online) Optimal |fc| as a function of the mass 
of the field for different vines of e — 6, 9, f 2, ... , 99 and p — 
10, 100, 1000 showing rapid saturation in e. For higher e this 
curves are completely insensitive to e variations, being very 
little sensitive for smaller values e < 10 




FIG. 7: 5g^''(e): Maximum entanglement achievable (optimal 
m and |fcl) as a function of e. It does not depend at all on p 



Although this method presents saturation when e ^^ oo 
(being most effective for e < 20) its insensitivity to p 
means that the optimal |fc| method gives us two indepen- 
dent methods for estimating p and e. In other words, all 
the information about the parameters of the expansion 
(both volume and rapidity) is encoded in the entangle- 
ment for the optimal frequency | fc | . 



C. Interpretation for the dependence of Se on |fc| 

We have seen (figure 1) that for bosons a monotonically 
decreasing entanglement is observed as |fc| increases. By 
contrast, in the fermionic case we see that there are priv- 
ileged |fc| for which entanglement creation is maximum. 
These modes are far more prone to entanglement than 
any others. 

To interpret this we can regard the optimal value of | fc | 
as being associated with a characteristic wavelength (pro- 
portional to |fe|~^) that is increasingly correlated with 
the characteristic length of the universe. As p increases 
the peak of the entanglement entropy shifts towards 
higher |fc|, with smaller characteristic lengths. Intu- 
itively, fermion modes with higher characteristic lengths 
are less sensitive to the underlying space-time because 
the exclusion principle impedes the excitation of 'very 
long' modes (those whose |fc| — >■ 0). 

What about small |fc| modes? As shown in [41] and 
in figure 1, for bosons the entanglement generation is 
higher when |fe| — )■ 0. This makes sense because modes 
of smaller | fc | are more easily excited as the space-time ex- 
pands (it is energetically much 'cheaper' to excite smaller 
|fc| modes). For fermions entanglement generation, some- 
what counterintuitively, decreases for |fc| -^ 0. However 
if we naively think of fermionic and bosonic excitations in 
a box we can appreciate the distinction. We can put an 
infinite number of bosons with the same quantum num- 
bers into the box. Conversely, we cannot put an infinite 
number of fermions in the box due to the Pauli exclusion 
principle. This 'degeneracy pressure' impedes those 'very 
long' modes (of small |fc|) from being entangled by the 
underlying structure of the space-time. 



V. CONCLUSIONS 

We have shown that the expansion of the universe 
(in a model 2-dimensional setting) generates in entangle- 
ment in quantum fields that is qualitatively different for 
fermions and bosons. This result is commensurate with 
previous studies demonstrating significant differences be- 
tween the entanglement of bosonic and fermionic fields 
[6, 7, 9, 16]. 

We find that the entanglement generated by the ex- 
pansion of the universe as a function of the frequency 
of the mode in the fermionic case peaks, while in the 
bosonic case it monotonically decreases. For bosons the 
most sensitive modes are those whose |fe| is close to zero. 
However for fermions modes of low |fc| are insensitive to 
the underlying metric. There is an optimal value of |fc| 
that is most prone to expansion-generated entanglement. 
This feature may be a consequence of the Pauli exclusion 
principle, though we have no definitive proof of this. 

We have also demonstrated that information about the 
spacetime expansion parameters is encoded in the en- 
tanglement between fermionic particle and antiparticle 
modes of opposite momenta. This can be extracted from 



the peaked behaviour of the entanglement shown in fig- 
ure 2, a feature absent in the bosonic case. Informa- 
tion about the rapidity of the expansion (p) is codified in 
the frequency of the maximahy entangled mode, whereas 
the information about the volume of the expansion (e) is 
codified in the amount of entanglement generated for this 
optimal mode. As e tends to infinity the maximum possi- 
ble S"^^^ in the optimal mode approaches the maximally 
entangled state. 

Hence the expansion parameters of spacetime are bet- 
ter estimated from cosmologically generated fermionic 
entanglement. Furthermore, these results show that 
fermionic entanglement is affected by the underlying 
spacetime structure in a very counterintuitive way and 
in a radically different manner than in the bosonic case. 
The manner and extent to which these results carry over 
to rf-dimcnsional spacetime remains a subject for future 



study. 
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